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Abstract. In this paper we consider the relation between the volume deceleration 
parameter obtained within the Buchert averaging scheme and the deceleration 
parameter derived from the supernova observation. This work was motivated by 
recent findings that showed that there are models which despite A = have volume 
deceleration parameter q""^ < 0. This opens the possibility that backreaction and 

^ ' averaging effects may be used as an interesting alternative explanation to the dark 

C^ . energy phenomenon. 

We have calculated q"°^ in some Lemaitre-Tolman models. For those models which 
CN ' are chosen to be realistic and which fit the supernova data, we find that g"°' > 0, 

^ i while those models which we have been able to find which exhibit 9"°' < turn out to 

be unrealistic. This indicates that care must be exercised in relating the deceleration 



l/^ . parameter to observations. 
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^ ■ 1. Introduction 



Accelerated expansion, modeled by a positive cosmological constant, is an essential 
element of the current standard cosmological model of the Universe. The accelerated 
expansion was originally motivated by supernova observations [T] and is supported 
by many other types of cosmological observations. Observational data is, in modern 
cosmology, analyzed almost exclusively within the framework of homogeneous and 
isotropic Friedmann models [2]. This analysis leads to the Concordance model, which 
provides a remarkably precise fit to cosmological observations. In this situation, if the 
Ehlers-Geren-Sachs theorem [3] and 'almost EGS theorem' |4] are invokeqj, then it 
seems that an assumption of large scale homogeneity of the Universe can be justified. 

I These theorems imply that if anisotropies in the cosmic microwave background radiation are small 
for all fundamental observers then the Universe is locally almost spatially homogeneous and isotropic. 
However, as shown in [5] the almost Robertson-Walker geometry also requires the smallness of the 
Weyl curvature. 
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This on the other hand implies that the Universe must be filled with dark energy which 
currently drives the acceleration of the Universe. 

However the Concordance model is not the only one which can fit cosmological 
observations. Anti-Copernican inhomogeneous models which assume the existence of a 
local Gpc scale void also fit cosmological observations [6] (see ||7j for a review). Moreover, 
on small and medium scales our Universe is not homogeneous. Therefore, one may ask 
whether Friedmann models can describe our Universe correctly. In particular, it is 
important to ask what is the best way to fit a homogeneous model to a realistic and 
inhomogeneous Universe. This problem, known as the fitting problem, was considered 
by Ellis and Stoeger [8]. In considering the fitting problem, it becomes apparent that 
a homogeneous model fitted to inhomogeneous data can evolve quite differently from 
the real Universe. The difference between evolution of homogeneous models and an 
inhomogeneous Universe is caused by backreaction effects, due to the nonlinearity of 
the Einstein equation. Unfortunately, in the standard approach, the backreaction is 
rarely taken into account - in most cases when modelling our Universe on a local scale 
Newtonian mechanics is employed and on large scales the Friedmann equations (or linear 
perturbations of Friedmann background) are used [9]. Such an approach to cosmology 
is often encouraged by the "no-go" theorem which states that the Universe can be very 
accurately described by the conformal Newtonian metric perturbed about a spatially 
flat background, even ii6p/p ^ 0. In such a case the backreaction is negligible [TOlfTT]. 
However, the results obtained by van Elst and Ellis [12] and recently by Kolb, Marra 
and Matarrese [I3] show that the application of "no-go" theorem is limited. Therefore, 
one should be aware that in the absence of an analysis of the backreaction and other 
effects caused by inhomogeneities in the universe, there remains the possibility that 
the observed accelerated expansion of the Universe is only apparent [Tl|. The direct 
study of the dynamical effects of inhomogeneities is difficult. Due to the nonlinearity 
of the Einstein equations, the solution of the Einstein equations for the homogeneous 
matter distribution leads in principle to a different description of the Universe than 
an average of a inhomogeneous solution to the exact Einstein equations (even though 
inhomogeneities when averaged over a sufficiently large scale might tend to be zero). 

Neither the analysis of the evolution of a general matter distribution nor the 
numerical evolution of cosmological models employing the full Einstein equations are 
available at the level of detail which would make them useful for this problem. There 
are currently several different approaches which attempt to take backreaction effects 
into account. One approach is based on exact solutions - see for example [15]. Another, 
and more popular approach is based on averaging. 

In the averaging approach to backreaction, one considers a solution to the Einstein 
equations for a general matter distribution and then an average of various observable 
quantities is taken. If a simple volume average is considered then such an attempt 
leads to the Buchert equations [16]. The Buchert equations are very similar to the 
Friedmann equations except for the backreaction term which is in general nonvanishing, 
if inhomogeneities are present. For a review on backreaction and the Buchert averaging 
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scheme the reader is referred to [I7l [I8]. Within this framework and using spherically 
symmetric inhomogeneous models Nambu and Tanimoto |19| . Paranjape and Singh |20| . 
Kai, Kozaki, Nakao, Nambu, and Yoo [2T], Chuang, Gu, and Hwang [22], provided 
explicit examples that one can obtain negative values of the volume deceleration 
parameter even if A = 0. Another interesting example was presented by Rasanen 
[TTl [23j where it was shown that the total volume deceleration parameter of two isolated 
and locally decelerating regions can also be negative. 

There are however important ambiguities in the application of an averaging 
procedure. The average itself not only depends on a choice of volume but also on a 
choice of time slicing. This is very crucial in cosmology. Once inhomogeneities are 
present the age of the Universe is not everywhere the same. Namely, the big bang in 
inhomogeneous models is not a single event, so the average taken over a hypersurface of 
constant cosmic time t is different from the average taken over a hypersurface of constant 
age of the Universe t — ts [24]. Moreover, the results of the averaging procedure vary if 
the discrepancy between the average cosmic time and the local time is introduced (the 
local time is the time which is measured by local clocks; the cosmic time is the time 
which appears in the averaged homogeneous model). This phenomenon was studied by 
Wiltshire [25], and has been used in an ambitious alternative concordance model. The 
model proposed by Wiltshire introduces some additional assumptions which allow to 
some extent a comparison of averaged quantities with observations. Such a comparison 
shows quite good agreement with observations, [26|- Thus, while serious fundamental 
questions remain concerning Wiltshire's approach, it is another example of an approach 
where one does not need dark energy to fit cosmological observations. 

The averaging procedure is also gauge-dependent. For example using different gauge 
one can obtain that the backreaction mimics not dark energy but dark matter [27|. The 
averaging schemes, therefore, in the literature have been criticized, and their inherent 
ambiguities (and in some cases obscurity) have been discussed, cf. e.g. [ID]. A key 
point is that it is far from obvious if the average quantities, such as the acceleration 
of the averaged universe are really the quantities which are measured in astronomical 
observations. In particular, an operational analysis is to a large extent lacking in the 
discussions of averaging. Thus, it is important to test the averaging procedures with 
the exact and inhomogeneous solutions of the Einstein equations. Within exact models 
each quantity can easily be calculated and then compared with its averaged counterpart. 
This paper aims to perform such an analysis within the Lemaitre-Tolman model. 

The structure of this paper is as follows. Buchert's averaging procedure is presented 
in section ^ and some background on the Lemaitre-Tolman model is given in section [3l 
The volume and distance deceleration parameters are introduced in section [H Finally, 
in section [5l we discuss the relation between the deceleration parameters, supernova 
observations and models of cosmic structures. 
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2. The Buchert scheme 

If the averaging procedure is applied to the Einstein equations, then for irrotational and 
pressureless matter the following equations are obtained [T6j 



3- = -AnG{p) + Q, (1) 

a 

3^ = 87^G(p)-i(7^)-^Q, (2) 

Q^l{{Q')-{er)-2{a'), (3) 

where (7^) is an average of the spacial Ricci scalar ^^^7^, 9 is the scalar of expansion, 
a is the shear scalar, and ( ) is the volume average over the hypersurface of constant 
time: {A) = (J cPx\/—h)^^ J (Px\/—hA. The scale factor a is defined as follows: 

a = iV/VoY/', (4) 

where Vq is an initial volume. 

Equations ^ and ([2]) are very similar to the Friedmann equations, where Q=0, 
and p and 7^ depend on time only. In fact, they are kinematically equivalent with a 
Friedmann model that has an additional scalar field source [28]. However the Buchert 
equation do not form a closed system. To close these equation one has to introduce 
some further assumptions [16]. As can be seen from (l3|) if the dispersion of expansion 
is large, Q can be large as well and one can get acceleration (a > 0) without employing 
the cosmological constant. 

3. The Lemaitre-Tolman model 



The Lemaitre-Tolman models [3^ is a spherical symmetric, pressure free and 
irrotational solution of the Einstein equations. Its metric is of the following form 

d.^ = cMt^ - ^1^ dr^ - R^it, r)d^^ (5) 

where dQ"^ = d6'^ + sin^^d0^. Because of the signature (+, — , — , — ), the E{r) function 
must obey E{r) > —1/2. Prime ' denotes dr. 

The Einstein equations reduce, in A = case, to the following two 

, , . 2M'(r) 

^^^^^'^' - R^ir,t)R'lr,ty ^'^ 



R{r,t) 



^-R\r,t) = 2E{r) + ^^, (7) 



§ The pressure free and irrotational solution of the Einstein equations for spherically symmetric space- 
time is often called the Tolman, Tolman-Bondi, or Lemaitre-Tolman-Bondi model. However, it is more 
justified to refer to this solution as to the Lemaitre-Tolman model (cf. |29j). 
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where M{r) is another arbitrary function and k = SttG/c'^. Dot ' denotes dt. 

When R' = and M' ^ 0, the density becomes infinite. This happens at shell 
crossings. This is an additional singularity to the Big Bang that occurs at -R = 0, M' ^ 0. 
By setting the initial conditions appropriately the shell crossing singularity can be 
avoided (see [STj for detail discussion). 

Equation (JT]) can be solved by simple integration: 

c\t-tB{r)\, (8) 



2E + ^ 



where t^ appears as an integration constant and is an arbitrary function of r. This 
means that the big bang is not a single event as in the Friedmann models, but occurs 
at different times at different distances from the origin. 
The scalar of the expansion is equal to 

^ R li , ^ 

The shear tensor is of the following form: 

c^% = ^§-f jdiag(0,2,-l,-l), (10) 

thus a^ = {l/2)a^i3a^f^ = {1/3){R'/R' - R/Rf. 

The spacial Ricci scalar in the Lemaitre-Tolman is equal to 

4. The apparent and average acceleration 

The deceleration parameter within the Friedmann models is defined as 

da 
? = --, 12 

where a is the scale factor. By analogy we can define the deceleration parameter which 
is based on the averaging scheme. Substituting (I4|) into f lT2] ) and using ([T]) and ([2]) we 
get 

vol _ -47rG'(p) + Q 

^ 87^G(p)-i(7^)-iQ• ^'-'' 

We refer to this deceleration parameter as the volume deceleration parameter, q™^ since 
it is positive when the second derivative of volume is negative and negative when the 
second derivative of volume is positive (and of sufficiently large value) . 

On the other hand one can introduce a deceleration parameter defined relative to 
the distance. Within homogeneous models the distance to a given redshift is larger 
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for accelerating models than for decelerating ones. Taylor expanding the luminosity 
distance in the Friedmann model we obtain 

1 d^Di 



D.. '^^ 



dz 



z + 



2 dz"^ 



z' + 0{z-') 



z=0 

= ^z + -^il-q)z' + Oiz'). (14) 

Ho /-Ho 

Employing a similar procedure in the case of the Lemaitre-Tolman model we get 
cR' cR'l^ R'R! cR" cR" , ., ^^, .,, 

R ^R' \ i?'2 R'R' R'2 I ^ ^' ^ ' 



Thus by comparing (iTSll with (IT411 . the Hubble and the deceleration parameter in the 
Lemaitre-Tolman model can be defined as 

rrdis _ R' ^dis _ R'R' cR" ^ cR" 

The above quantities are defined at the origin (r = 0). However, following the Partovi 
and Mashhoon [32] we can extend the above quantities to any r. Then, the coefficients 
of Taylor expansion are 

^= ^R + R^ + R'^ 

dz dz dz 

d^Dr ■ dt ,dr ■■ f dt\ -.dtdr „ f dr\ 

-2R + AR— + AR!— + R\ — ] + 2R'—— + R" { — ] 



and we obtain 



dz'^ dz dz \dz J dz dz \dz J 

;^dH d^r 

dz"^ dz 



+ R—.+R'—2^ (17) 



H"^-^ = (^] ' , q''' = 1 - H''^^. (18) 

\ dz J dz^ 

We refer to this deceleration parameter as the distance deceleration parameter. Although 
of physical importance is the luminosity distance and its ability of fitting the supernova 
data, the g"^** is of great usefulness. It allows us, without solving the geodesic equations, 
to easily check whether a considered model can be use to fit supernova data. As we will 
see in the next section, models which fit supernova data have at least in some regions 

5. Connection between deceleration parameter and observations 

Let us first focus on supernova observations. There is already a considerable literature 
on inhomogeneous models which are able to fit the supernova observations without the 
cosmological constant [6j. We shall examine four such models in this section. For each 
of these models we shall calculate the volume and distance deceleration parameters and 
compare with each other. The four models to be considered present a very good fit 
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Figure 1. The Residual Hubble diagram for models 1-4. The black dashed line 
presents Am for the ACDM model. 



to supernova data. The supernova data consists of 182 supernovae from the Riess gold 
sample [3_3J. The x^ test for models 1-4 is respectively 183.6, 184.3, 164.7, and 178.5 (for 
comparison the x^ of fitting the ACDM model is 165.3). The residual Hubble diagram 
for these models is presented in figured The deceleration parameters for models 1-4 are 
presented in figure [2l Left panel presents the distance deceleration parameter [as defined 
by ( fTSl) - where dt/dz and dr/dz were calculated for the radial geodesic]. The distance 
deceleration parameter is positive at the origin, but soon becomes negative. Moreover, a 
very similar shape is obtained if instead g'^*'' [as defined by (ITSl) ] gg*"* [^^ defined by (IT6l) ] 
is used. Thus, g'^*'^ (or even ^q**, if treated as a function of r) can be regarded as a useful 
test to check if a given model is able to fit supernova data. However, the most significant 
is that the volume deceleration parameter which is presented in the right panel of figure 
[2] is strictly positive. Thus, the ability of reproducing the supernova data does not 
require that the volume deceleration parameter is negative. This raises the question 
whether the average acceleration has any relation with the observed acceleration of the 
Universe; and if yes, are models with average acceleration also able to fit supernova 
data? 

Let us now focus on models of cosmic structures. It was recently shown that 
using a perturbative approach, backreaction cannot explain the apparent acceleration 
[34] . However, because of large density fluctuations within cosmic structures, results 
obtained in terms of the perturbation framework might be questionable. Moreover, 
in view of the fact that there are known examples of exact inhomogeneous models 
with negative volume deceleration parameter and A = 0, it is worthwhile to check if 
realistically evolving models of cosmic structures can have negative values of deceleration 
parameter. First, let us consider a model of galaxy clusters with the Navarro-Frenk- 
White density distribution [35] (left panel of flgure [3|). Although, the NFW proflle 
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Figure 2. The distance deceleration parameter (left panel) and the volume 
deceleration parameter (right panel) for models 1-4. 



0.9 

0.7 

0.5 

0.3 

0.1 



-0.1 



5a, 



5b 



0.2 0.4 0.6 0.8 

R [Mpc] 



0.2 0.4 0.6 0.8 

R [Mpc] 



Figure 3. The current density distribution (left panel) and deceleration parameter 
(right panel) for model 5. 



describes virialized system^ the use of this profile will prove to be very instructive. 
The average deceleration parameter g^°' for model 5 is presented in the right panel of 
figure [3l As can be seen in this case the deceleration parameter is positive (curve 5a). 
However, it is possible to modify this model so that the q'"°^ becomes negative - curve 5b 
in the right panel of figure [3l This was obtained by choosing the E function which is of 
large positive value (for details see Appendix). However, after such a modification this 
model becomes unrealistic. Specifically, the age of the Universe in this model becomes 
unrealistically small. The bang time function t^ in this model is of large amplitude, 
around 11.44 x 10^ y. This means that the actual age of the Universe in this model is 
approximately a few hundreds of thousand years. 

II The Lemaitre-Tolman model which evolve from smooth density profile at last scattering to a high 
value profile like the NFW profile is always characterized by a collapse - central region within this 
models are at the current instant collapsing. Thus such systems cannot be considered as virialized 
systems. 
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Figure 4. The current density distribution (left panel) and volume deceleration 
parameter (right panel) for models of cosmic structures (models 6, 7). 
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Figure 5. The current density distribution (left panel) and deceleration parameter 
(right panel) for model 8. 



Now let us examine the volume deceleration parameter within models of cosmic 
voids and superclusters. Figure [4] presents density distribution of realistically evolving 
cosmic structures (void - curve 6, supercluster - curve 7). It can be seen from the 
right panel of figure [4] that the volume deceleration parameter within these models is 
positive. As above, we can modify our models in such a way that the volume deceleration 
parameter is negative, but again this leads to a very large amplitude of t^. For example, 
in model 8 whose density and the volume deceleration parameter are presented in figure 
[Qlf l the volume deceleration parameter is negative. However, the bang time function in 
model 8 is of amplitude ^ 11 x 10^ y, which leads to unrealistically small age of the 
Universe. 

% Employing a model of qualitatively similar features as model 8, Hossain [36j showed that the observer 
situated at the origin in order to successfully employ the Friedmann model has to assume the existence 
of dark energy. 
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6. Conclusions 

In this paper we have studied the relation between the volume deceleration parameter 
obtained within the Buchert averaging scheme and the deceleration parameter derived 
from the observations of supernovae. This work was motivated by recent results showing 
there there are models which despite A = and average expansion rate is accelerating, 
i.e. a > [where a is defined by relation (I!])]. This opens the possibility that 
backreaction and averaging effects may be used as an interesting alternative explanation 
to the dark energy phenomenon. 

We have compared the quantities obtained within the exact and inhomogeneous 
models with their average counterparts. We focused on the supernova observations 
and models of cosmic structures. For this purpose the Lemaitre-Tolman model was 
employed. It was showed numerically that the averaging of models which fit the 
supernova observations does not lead to volume acceleration (a < for these averaged 
models and hence g"°' > 0). It was also shown that realistically evolving models of 
cosmic structures have also g^°' > 0. It was possible to modify these model in such 
a way that after the averaging q^"'' < 0. This was obtained by choosing E function of 
positive amplitude - as was recently proved by Sussman [37] this is a necessary condition 
to obtain q^°^ < 0. However, in models with realistic density distribution, in such a cases, 
£" > 1 > M/R ^ 10"'^ - 10"^, hence as seen from ^ ts ^ t (to remind c x 10^°y ^ 3 
Gpc). Thus, within such models the age of the Universe is unrealistically small. 

Our analysis has been performed in the limited class of Lemaitre-Tolman models, 
which due to their spherical symmetry are arguably too simple to give a full 
understanding of averaging and backreaction problems. However, within this class, 
we conclude that the volume deceleration parameter q^°'- is not a quantity which can be 
directly related to observations. 

It is possible that the volume deceleration parameter g""' becomes negative only 
after averaging over the scales which are larger than 100 Mpc. On such large scales the 
structure of the Universe becomes too complicated to be fully described by spherically 
symmetric models. However, it is intriguing that models which fit the supernova 
observations and for which the distance deceleration parameter, g'^*'^, is negative have 
still g^"' > 0. This suggest that the volume deceleration q^°^ does not have a clear 
interpretation in terms of observable quantities. It does not, of course, mean that 
averaging and backreaction effects cannot potentially be employed to explain the 
phenomenon of dark energy. However, our work here indicates that such a potential 
solution of the dark energy problem should be based upon different methods than 
those related to volume deceleration parameter. Rather than showing that q'"°^ < the 
averaging approach should explain observations - reproduce correct values of distance 
to supernovae, correct shape of the CMB power spectrum, etc. An interesting, quasi- 
Friedmannian approach, was recently suggested in [38]. In this approach backreaction 
is modeled in terms of the morphon field [28]. In such a case a Universe is describe by a 
homogeneous model with the spatial curvature being just a function of time. As shown 
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in [38] such approach lead to an agreement with supernova and CMB data without the 
need for dark energy, but requires q^°^ < 0. 
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Appendix A. Model specification 

There are three arbitrary functions of the radial coordinate in the Lemaitre-Tolman. 
However only two functions are independent and the third one is specified by the choice 
of the radial coordinate. Models considered in this paper are defined as follows: 

(i) Model 1 and 2 

The radial coordinate is chosen as the present day value of the areal distance 
r := -Ro- Models 1 and 2 are specified by the present day density distribution 
and the time bang function. The density distribution is parametrized by 

,2 



p{to, r) = pb 



r 
1 + 5p- 5pexp ( 



a2 



(A.l) 



where p^ = fi„ x {3H^)/{87tG), Qm = 0.27, Hq = 70 km s'^ Mpc'^ In model 1 
ps = 1.9, cr = 0.9 Mpc, and in model 2 ps = 1.5, a = 0.5 Mpc. In these models 
the big bang is assumed to occur simultaneously at every point, i.e. ts = 0. The 
functions M and E are then calculated and using eqs. ((61) and ([8]) respectively. 
The time instants as well as background density pb in all models (1-8) is chosen as 
density of a Friedmann model {flm = 0.27, Hq = 70 km s~^) and time instants are 
calculated using the following formula [9j: 

oo 

ti^) = ^ [ "^"^ (A.2) 

z 

where Qk = 1 — ^m- The last scattering instant (tLs) is set to take place when 
z = 1089 and the current instant (to) when z = — tis = 4.98 x 10^ y, and 
to = 11.4421 X 10^ y. 

(ii) Model 3 and 4 
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As above the radial coordinate is chosen as the present day value of the areal 
distance r := Rq. These two models are defined by the current expansion rate, and 
assumption that p{to, r) = ph. The expansion rate is parametrized using 



HT{to,r) = - = Ho 



1-6h + SHexp [ 



(A.3) 



where E^bu = 9.6 km s"^ Mpc"\ a = 0.6 Mpc, and Ho6h = 12 km s"^ Mpc'^ 
0" = 1.2 Mpc for model 3 and 4 respectively. In these models density is assumed 
to be homogeneous at the current epoch. The function M is then calculated using 
the above relation and eq. (jTl). It should be noted that the Ht is one of several 
generalization of the Hubble constant, which in the Friedmann model is Hq = a/a. 
Apart from the transverse Hubble parameter, Ht, one can also define the radial 
Hubble parameter, Hji, [see eq. (IT6ll]. and the volume Hubble parameter defined 
as Hv = (1/3)0 = Hr + 2Ht. 

(iii) Model 5a 

The radial coordinate is chosen as the present day value of the areal distance, 
i.e. r := Rq. The model is defined by density distributions given at the present 
instant and at last scattering. The density distribution at the current instant is 
parametrized by 

p{to, r) = pb — -— , (A.4) 

{r/rs){l + r/rsy 

where 6 = 28170 and r^ = 191kpc. This is a Navarro, Frenk, and White galaxy 
cluster profile [35]. As can be seen this profile is singular at the origin but this 
problem can be overcome by matching the NFW profile with a singular-free profile 
as /(r) = —ar^ + b. 

The density profile at last scattering is assumed to be homogeneous, thus the areal 
distance at last scattering is: 

RitLs,r)=( -) . (A.5) 

The function M{r) is then calculated from eq. ([6|). Function E can be calculated 
by subtracting solutions of ([8|) for t^s and to (for details see [39]). The function E 
is presented in the left panel of figure lAli 

(iv) Model 5b 

The radial coordinate is chosen as the present day value of the areal distance, 
r := Rq. The model is defined by density distribution given by (IA.4I) and E of the 
following form 

E{r) = 10^ sin (lO-^Mpc"^) . (A.6) 



This profile is presented in the left panel of figure lAll and the bang time function 
tB in the right panel. 
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Figure Al. Left panel presents the function E{r) for models 5a and 5b. Please note 
that the y-scale in the upper part of the left panel is different than in the lower part. 
Right panel bang time function for model 5b. 



VI 



Models 6 and 7 

The radial coordinate is chosen as the value of the areal distance at last scattering 
instant, r := Rls- Model 6 and 7 are defined by the assumption that t^ = and 
the density distribution, which at last scattering is of the following form 



P(^L5, r)= py,\ 1 - (5 exp {-al r ) + 7 exp 



Ir 



d 



(A.7) 



where i = 1/kpc; 6 = 1.2 x 10^^ and 2 x 10"'^ for model 6 and 7 respectively; 
7 = 14.62 X 10"'^ and 8.03 x 10""^ for model 6 and 7 respectively; a = 0.01 and 
0.04 for model 6 and 7 respectively; c = 18 and 12 for model 6 and 7 respectively; 
and d = 6 and 5 for model 6 and 7 respectively. The bang time function for both 
these models is ts = 0. The mass function, M{r) is calculated from eq. ([6]), and 
the function E{r) is calculated from eq. ([8]). 

Model 8 

The radial coordinate is chosen as a present day value of the areal distance: r := Rq. 

Density distribution is of the following form 



p(to,r) = 6.2pfeexp(-4 x 10"*(£r)^) 



and the function E is 



Fir) 



—r) exp(10-3£r), 



(A.8) 



(A.9) 



which except for [exp (10~'^£r)] is the same as E{r) profile in the empty Universe. 
This profile is presented in the left panel of figure IA2I and the bang time function 
ts in the right panel. 
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